THE LEADING BEHAVIOR OF THE GROUND- STATE ENERGY 
OF HEAVY IONS ACCORDING TO BROWN AND RAVENHALL 



XIAO LIU 



Abstract. In this article we prove the absence of relativistic effects in leading 
order for the ground-state energy according to Brown-Ravenhall operator. We 
obtain this asymptotic result for negative ions and for systems with the number 
of electrons proportional to the nuclear charge. In the case of neutral atoms 
the analogous result was obtained earlier by Cassanas and Siedentop [2]. 
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1. Introduction 

Cassanas and Siedentop [4] have shown that the ground state energy of heavy 
atoms for the relativistic Hamiltonian of Brown and Ravenhall is, to leading order, 
given by the non-relativistic Thomas-Fermi energy. The relativistic Hamiltonian 
of Brown and Ravenhall is derived from quantum electrodynamics yielding energy 
levels correctly up to order a 2 Ry [3]. However, only the case as N — Z is considered 
in [4] , where N is the number of electrons and Z is the nuclear charge. This thesis 
describes the other two cases: N > Z and N/Z — X(constant). 



2. Definition of The Model 



Brown and Ravenhall [5] describe two relativistic electrons interacting with an 
external potential. The model has an obvious generalization to the iV-electron case. 
First we define 



(1) 



D r Z ■= a 



fp-Z\-\~ 1 



is the Dirac operator of an electron in the field of a nucleus of charge Z . Note that 
we are using atomic units in this paper, i.e., m e = H = e = 1. As usual, the four 
matrices ai, «2,Q?3 and /3 are the four Dirac matrices in standard representation, 
explicitly 

<T 0. 



<r denoting the three Pauli matrices 



CTi = 



1 

1 



ct 2 = 



-i 

1 



0-3 



1 

-1 



and 



(\ 

10 

0-10 

\0 -1/ 



Let Qat := AiLi(# 1/2 (^ 3 ) ® C 4 ) n f) N where 
(2) iojv := A 

the underlying one-particle Hilbert space is 



(3) 

where 
(4) 

for any ip € L 2 



S:= [x ( o,oo)Pc,o)](i 2 (K 3 )®C 4 ) 



X(o,oo) (D c ,o) 
C 4 , and 



1, D c ,oV) > 0; 
0, (V»,Dc,o^)<0 



i 
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Now we define the model as 
£ : £} N -> R 

(K\ N 

[) ^^(v>,(2p c , -c 2 -z/|x w |) 1/ + X! i^-x.r 1 )^). 

i/=l l<n<v<N 

As an immediate consequence of the work of Evans et al. [5] quadratic form £ 
is bounded from below, in fact it is positive (Tix [33 [36]), if k := Z/c < K C rit := 
2/ (tt/2+2/tt). According to Fricdrichs this allows us to define a self-adjoint operator 
B c ,N,z whose ground state energy 

(6) E(c, N, Z) := inf a(B c , N , z ) = M{£{^ € Q Nl = 1} 

is of concern to us in this paper. In fact the main result of this article is 
Theorem 1. When N > Z or N/Z — X(constant) , we have 
E(Z/n, N, Z) = E TF (N, Z) + o(Z 7/3 ) 
where E TF (N, Z) := inf{£ TF (p) \p > 0, / R3 p < N, p G L 5 / 3 (M 3 )}, 

dx + D{p, p) 



£tf(p) 



^7TFP(x) 5/3 - n p(x) 

5 x 



is the Thomas-Fermi functional, and n := Z/c < « cr it := 2/tt. 

In the following, we will assume that the ratio n € [0, ft cr it) is fixed. Note that 
according to [531 Formula (9c)] the minimal energy E^y{N,Z) fulfills the scaling 
relation 

(7) E TF {N, Z) = E TF (N/Z, 1)Z 7/3 . 

The article is structured as follows: first we show how the treatment of the 
Brown-Ravenhall model can be reduced from Dirac spinors (4-spinors) to Pauli 
spinors (2-spinors) . Then we separately give the upper and lower bounds for the case 
N > Z in Section 3 and for the case where N/Z is fixed in Section 4. Throughout 
the paper, we use the letter k for any constants independent of c, N, R, or Z. 

We now indicate, how to reduce to Pauli spinors. To this end we parameterize 
the allowed states: Any ip £ f) can be written as 

(8) W > 



a 



\ N c (p) 

for some u E \) := L 2 (R 3 ) <£> C 2 . Here, <x are the three Pauli matrices, 

p := -iV, E c (p) := (c 2 p 2 + c 4 ) 1 / 2 , N c (p) := [2E c (p)(E c (p) + c 2 )] 1 ' 2 . 
In fact, the map 

$ : fj -> Sj 



(9) f<S> lU 

$2" 



/ ge(p) + C 2 \ 
. = W £ (p) " 



embeds f) unitarily into ^ and its restriction onto i? 1 (E 3 ) (g> C 2 is also a unitary 
mapping to ^ n H 1 {R 3 ) ® C 4 (Evans et al. 0). 
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It suffices to study the energy as function of u 

N 

(10) £ o(®^ =1 $) : f\ fi^R. 

The one-particle Brown- Ravenhall operator B 1 for an electron the external elec- 
tric potential of a point nucleus acting on Pauli spinors is then (see Appendix [A| 

(11) B CtZ = E c (p) - Zipx - Z V2 
where 

(12) B CtZ := X{a,oo){Dcfi)D c , z 
and we have split the potential into 

(13) <pi := $;| • l" 1 ^, ^ 2 :=$5| • |- 1 $ 2 . 

As we will see the first part (pi is contributing to the nonrelativistic limit whereas 
the second part turns out to give the energy contribution that does not even affect 
the first correction term. 

3. Case I: N > Z 
In this section, we prove Theorem [T] for negatively charged atoms. 

3.1. Coherent States. We obtain the upper bound by constructing a trial density 
matrix in the Hartree-Fock functional for the Brown- Ravenhall operator. To this 
end we introduce spinor valued coherent states. 

Given functions /, / € ff 3 / 2 (R 3 ) and an element a — (p, q, r) of the phase space 
T := K 3 xR 3 x {1, 2}, we define coherent states in f) as 

(14) F a (x) := ((y5p : q ® e T )(x) := f(x - q) exp(ip • x)S T , a , 

where x — (x, a) £ M 3 x {1, 2} and the vectors e r are the canonical basis vectors 
in C 2 (see Lieb |2T] and Evans et al. [S]). We also define 

(15) ^(x) := (2^)- 3 / 2 %), 

where k S N. We will pick / depending on a dilation parameter. More specifically, 
we will choose 

(16) /(x) := 5R (x) := i?- 3 / 2 .g(i?- 1 x) 
(see [H Formula (11)]), 

'(2TTR)- 1 / 2 \\x\-R T 1 sin(7r | |x| - R \/R), R < |x| < 2R; 
0, otherwise, 



(17) /» := 



where R := Z~ s with 5 € (1/3, 2/3), R G K + and g € H 3 / 2 , spherically symmetric, 
normalized, and with support in the unit ball. 

The natural measure on T counting the number of electrons per phase space 
volume in the spirit of Planck is J r dil(a) := (2tt)~ 3 J dp Jdq^r=i- The essential 
properties needed are the following. For L 1 (r,dfi), let 

7i := J^Q(a)A(a)\F a )(F a \, 
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f r an(a)A(a) 



K+N-Z 
k=K+l 

K+Z 

72— ^ 2 I0jv-z+sWjv-z+sI> w/iere e fl := 
fe=i<"+l 

Then 

7 := 7i + 72 + 72 

and 

(18) 7l := J^dQ(a)A(a)\F a ){F a \ 

(see [H Formula (12)]) are trace class operators. 
We will pick 

(!9) A (a) X{(p,q)eR 6 | P 2 /2-v z (q)<o}(p,q) 

(see [H Formula (16)]), 

( 20 ) ^( a ) := X{(p,q)eR 6 |p 2 /2-V z (q)<0, |q|<A— J?}(P' ^ ' 

where Vz '■= Z/\ ■ \ — \ ■ \~ 1 * , here p^ Z ^ is the unique minimizer of the 

Thomas-Fermi functional 



(21) £ TF (p) := f 



^ 7 TFP(x) 5/3 - ||rp(x) 



dx + D(p, p) 



in the set of functions p € {p |p > 0, J R3 p < N, p e L 5/3 (R 3 )}. Here 
D(p,p):=\l I p(x)|x-y|- 1 p(y)dxdy 

^ ./R 3 ./R 3 

is the Coulomb scalar product. For fermions with g spin states per particle, 
7tf := (67r 2 /g) 2/3 ft 2 /(2m)(see Lieb [2TJ Formula (2.6)]), i.e., in our units, 7tf = 
(3^ 2 ) 2 / 3 /2. 

According to [25], E TF (N,Z) := inf£ TF (p) is strictly monotone decreasing for 
N < Z and constant for N > Z. Thus for N > Z, the minimizer p of the Thomas- 
Fermi functional (l2~Tj) coincides the minimizer in the case N — Z, which we denote 
by p z - Hence V z = Z/\ ■ \ - \ ■ \^ * p z for N > Z. 

By Lemma [T8l 

(22) < A < 1 =>■ < 7i < 1, 
according to the definition of A and the support sets of <pk, 

(23) < i < 1 < 7 < 1; 
by Lemma [T9l 



(24) tr 7 i = y JO(a) J 4(a) < N, 

(25) tr 7 = tr f! + tr 7 2 + tr 7~ 2 = iV. 
Using $ we can lift 7 to an operator on S) 

(26) 7$ := $7$*. 
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Note that 

(27) V^^^Z^V^iZ^q) 

(see Gombas jH] and [23]). Note also that ([57]) does not depend on N for TV > Z. 

3.2. Upper Bound on £^ F (j). We begin by noting that the Hartree-Fock func- 
tional - with or without exchange energy - bounds E(c, N, Z) from above. To be 
exact we introduce the set of density matrices 

(28) S aN := { 7 € e 1 ^) \E c (ph G S 1 ^), < 7 < 1, tr 7 = N} 
where G 1 ^) denotes the trace class operators on \). 

^hf : S 9N -> K 

(29) n 

7 M- tr[(£ c (p) - c 2 - Z/|x|) 7 »] + £>(p 7s ,p 7 J 

where - as usual - p 7 (x) := 7(x, x) is the density associated to 7 and -D is the 
Coulomb scalar product. By the analogue of Lieb's result [22] Formula (9)] and 
[20] (see also Bach p]) - which trivially transcribes from the Schrodinger setting to 
the present one - we have for all 7 £ Sqn 

(30) E(c,N,Z) <£« F ( 7 ). 
3.2.1. Kinetic Energy 0/71. 

Lemma 1. 77ie kinetic energy 0/71 does raoi exceed the kinetic energy 0/71, i.e., 

(31) tr[(F c (p) - c 2 )7 1$ ] < tr[(£ c (p) - c 2 ) 7l$ ]. 
Proof. According to 

(32) E c (p) - c 2 > 0, 
we have 

(33) 



tr[(£ c (p)-c 2 )(7 1$ - 7l$ )] = / / an(a)[A{a)-A{a)}[(E c {p)-c 2 )F a (x)}F a (x)dx 



r 



-lj(2n)- 3 J ^dpdq[(^ c (p)-c 2 ) 1 /2 Fa (x)][(i; c ( P )-c2)V2F a (x)]d^ 

p 2 /2-V z (q)<0 T = 1 
q|>-R--R 

= -(2tt)- 3 /" y'|(£; c (p)-c 2 ) 1 / 2 J F 1 Q (^| 2 dxdpdq<0. 



p' I /2-V ? (q)<0 
q|>H-R 



Thus tr[(£7 c (p) - c 2 )7i $ ] < tr[(£ c (p) - c 2 ) 7l$ ]. □ 
3.2.2. External Potential of 71. 

Lemma 2. For any e > 0, £/iere exists R large enough such that 

tr[(— Z/|x|)ti $ ] < tr[(— Z/|x|)7i $ ] + e , obviously for any fixed R, we have 

tr[(-Z/|x|)f 1$ ] > tr[(-Z/|x|) 7l$ ], i.e., 

(34) tr[(-Z/|x|)7i*] — > tr[(-Z/|x|) 7 i*], £ -> 00. 
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Proof. According to Lieb [2lJ Formula (2.18)], we can write the Thomas-Fermi 
equation as 

(35) 7tfPz /3 (x) = [V z ( X )-u'} + , 

where it' > is some constant and for t £ R, we set [t]+ := max{t, 0}. Then for 
N > Z, the Thomas-Fermi potential Vz '■= Z/\ • | — pz * | • | _1 can be written as 

(36) 7tfP z /3 = V z 
(see, e.g., Gombas [H]). 

Using Kato's inequality (see [HI Formula (2.9)]), passing from £ to 4 + p and 
taking into account that 

|ff(£)l 2 0£l 2 +2|£l + i] 

is integrable because g € _ff 3 / 2 (R 3 ), we can get 

(37) tr[(-Z/|x|)(7U -7i $ )] < \ Z J J dQ(a)[A(a) A(a)]\F a (^\m 

= ^zj (2rr)- 3 J ]Tdpdq|F Q (0| 2 |£|d£ 

p 2 /2-V z (q)<0 T_1 
|q|>S-fl 

<kZ J dSl(a)[A(a) - A{a)]R~ A J \g{£)\ 2 \£ + Rp\d£ 
<kZ J an(a)[A(a)-A(a)}R- A (l + R\p\) 
<kZ 1+4S J dpdq(l + i?|p|) 

p 2 /2 _ 7TFp 2/3 (q) < 
|q|>H-H 

< kZ 1+iS J [p z (q) + i? (0 2 /3 (q)]dq — ► 0, when R^oo. 

\q\>R-R 

The last step is according to absolute continuity of Lebesgue integral. 

Thus tr[(-Z/|x|)fi # ] — > tr[(-Z/|x|) 7l$ ], when R -> oo. □ 

3.2.3. T/ie Electron- Electron Interaction 0/71. 

Lemma 3. 77ie electron- electron interaction of 71 does noi exceed the electron- 
electron interaction 0/71, i.e., 

( 38 ) ^C/W^yiJ ^ D (P-yi*>P-n*)- 
Proof. 



(39) D(p^ lt ,Pf 1# ) 
1 

~ 2 
1 



< 
- 2 



dft(a)i(a)|<FF Q (x)| 2 |a; - y^ 1 y dfi(a)i(a)|$ J F Q (y)| 2 dxdy 
dfi(a)A(a)|$F Q (x)| 2 | a; - ^ dr>(a)A(cv)|$F Q (y)| 2 dxdy 
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□ 

3.2.4. Kinetic Energy 0/72. We introduce the set of density matrices 

(40) S N := {71 G e\t)) |£ c (p)7i G ©'(fl), < 71 < 1, tr 7 i < iV} 
where S 1 (f)) denotes the trace class operators on \). 

Lemma 4. For any e > 0, there exists a large enough K such that 
tr[(£ c (p) - c 2 ) 72<1 ,] < e , i.e., 

(41) tr[(£ c (p) - c 2 ) 72$ ] — ► 0, A" -> 00. 
Proof. By concavity we have 

(42) £ c (p)-c 2 <± P 2 = -±A, 

which implies that the Brown-Ravenhall kinetic energy is bounded by the non- 
relativistic one, i.e., for all 72 € SW with — A72 € © 1 (f)) 

(43) tr[(£ c (p) - C 2 ) 72 ] <tr(-±A 72 ). 
Then 

1 K+N-Z 

(44) tr[(S c (p)- C 2 ) 72$ ]<tr(-iA72)-2 E ( 2 Vl|V/l| 2 

1 00 1 fl^+i 

<^ E 4- fc ||V/|| 2 = i^^||V/f 

k=K+l 4 

= |(J) K+1 ||V/|| 2 ^0, tuten tf^oo. 

Thus tr[(£ c (p) - c 2 )7 2<1 ,] — > 0, when if -)• 00. □ 

3.2.5. External Potential 0/72- 

Lemma 5. The external potential 0/72 does not exceed zero, i.e., 

(45) tr[(-Z/|x|) 72 *] < 0. 
Proof. 

(46) 

K+N-Z K+N-Z 

tr[(Z/|x|) 72 *] = Z = Z E J l $ ^( x )| 2 O dx ^ °- 

□ 

3.2.6. The Electron-Electron Interaction 0/72. 

Lemma 6. For any e > 0, there exists a large enough K such that 
(47) D{p J2 ^ p 72 J — > 0, K ^ co. 
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Proof. According to [2 Equation (12)], we can get 
(48) £>(p 72 *,p 72 J < Jtr| 72$ |tr(^A| 72$ |) 

K+N-Z 

<\{n-z) (*hAP\*te) 

k=K+l 
K+N-Z 

<\{N-Z) I l^(P)l 2 |Pl d P 

~k=K+l 

oo „ 

k=K+l 

= 1(N - Z)(~) K J \f(Z)\ 2 \m ^ 0, when K ^ oo. 



Thus D(p 72# ,p 72 J — > 0, when K -> oo. □ 

3.2.7. The Total Energy of£§ F ( , y). We define the reduced Hartree-Fock functional 
on 71 as following 

: S N -> E 

(49) 

7i h+ tr[(E c (p) - c 2 - Z/|x|) 7 i*] + D(p 71# ,p 71 J 

where - as usual - p 7l is the density associated to 71 and D is the Coulomb scalar 
product. Gathering our above estimates allows us to get Theorem [2] 

Theorem 2. The reduced Hartree-Fock functional ofj does not exceed the reduced 
Hartree-Fock functional 0/71, i.e., 

(50) £&( 7 ) < £hV(7i), 
w/ien i? cm<i K are tending to infinity. 

Proof. According to the definition of 72 and Appendix [C] we get that as R tends 
to infinity, tends to zero, and thus £^(72) tends to zero. Using the results 
obtained in Lemmata [T] through [6l we get for any e > 0, there exist large enough 
R and K such that 

(51) £* F ( 7 ) = tr[(£ c (p) ~c 2 ~ Z/|x|) 7 «] + £>(p 7s ,p 7 J 

< tr[(£ e (p) - c 2 - Z/|x|) 7 i # ] + D(p 71# , p 7l J + e = £*,( 7l ) + e. 

□ 

3.3. Upper Bound. We begin by noting that the Hartree-Fock functional - with 
or without exchange energy - bounds E(c, N, Z) from above. To be exact we 
introduce the set of density matrices (|40p . where 6 1 (f)) denotes the trace class 
operators on f). We define the reduced Hartree-Fock functional of 71 as (|4l?|) . By 
Theorem [5] and the analogon of Lieb's result [HI [30] (see also Bach pQ) - which 
trivially transcribes from the Schrodinger setting to the present one - we have for 
all 71 e Sjv 

(52) S( C ,7V,Z)<^ F ( 7 )<f« F ( 7l ). 
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3.3.1. Kinetic Energy. (|42|) implies that the Brown-Ravenhall kinetic energy is 
bounded by the non-relativistic one, i.e., for all 71 6 Sn with — A71 € © 1 (ll) 

(53) tr[(£ c (p) - c 2 ) 7l ] < tr(-±A 7l ). 

Insertion of 71 (see Equations (fl"4|) . (fTo) . (H"8l) . and (IT9)) ') turns the right hand side 
into the Thomas-Fermi kinetic energy modulo the positive error (tr O'l)!! V<?|| 2 i? — 2 
(see Lieb [HI Formula (5.9)]). 
In fact, we choose 

/(x) = (27ri?)- 1 / 2 |x|- 1 sin(7r|x|/i?) 

(see Lieb [HI Formula (5.11)]). Because of /(x) = i?- 3 / 2 .g(i?- 1 x), we know 

.g(x) = (2^)- 1 /2| xr i sin( ^| x | ) . 

Let 77 = i? _1 x, we can calculate 

V/ = i?- 3 / 2 V.g(i?- 1 x) = R-^ Vvg{v) ^L = i? -3/2 v ( } 1 x 

ox it |x| 

|V/| =J R- 5 / 2 |V^(r,)|, | |V/| 2 d 3 x = I i?- 5 |V r)5 ( J 7)| 2 (i? 3 d 3 J7 )=i?- 2 ||V 3 || 2 . 
Thus we obtain 

(54) trp c (p) - c 2 ) 7 i] < jj 7 TF J 4 /3 (x)dx + Zi?- 2 ||V.g|| 2 . 

3.3.2. External Potential. Since — Ztr(</3 2 7i) is negative, we can and will estimate 
this term by zero. This estimate will be good, if this term is of smaller order. 
Although, logically unnecessary for the upper bound, it is interesting to see that 
if2 does indeed not significantly contribute to the energy, if 71 is chosen as above. 
Moreover, the proof will be also useful for the proof of Lemma [8] 

Lemma 7. For our choice 0/71 = J r dfl(a)A(a)\F a )(F a \ and S e (1/3,2/3) we 

have 

(55) 



J J ^ i^-i'i 2 A^v c (0 ( 



< Ztr(^27i) < kZ J dCl(a)A(a) 
Proof. We begin by estimating the expectation of ip2 on the coherent state ((T 



(56) 0< (F a ,ip 2 F a ) <k jj d£d£ 
< kc~ 2 R- 3 JJ dgdg' ^ffl \t+Rp\\?+Rp\ < kc- 2 R- 3 (l+R\p\+R 2 \p\ 2 ). 
Here, we have used that N c (£) > ^/2c 2 and, at the last step, that 

»fp^iiri + i^i + iri + i] 
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is integrable in £ and £' because g € H 3 / 2 (R 3 ). Thus according to (|56| . we get 



(57) < Zti{ m ) = Z I dn(a)A(a)(F a ,ip 2 F a ) 

^ k ^E J an(a)A(a)(l + R\p\+R 2 \p\ 2 ) 



< k 



c 2 R 3 



Z + R / dq 



R 2 / dq 



Z 4/ 3yi(z l/3 q) 



5/2 



Z 4/3 Fi(z l/3 q) 

= 0(Z 3S + z 2 ^ 3+2S + z 4 /' 3+s ) 

(see [1 Formula (27)]). □ 
Lemma 8. For our choice 0/71 and 5 G (1/3,2/3) we /love 
(58) |^tr[(| - j- 1 - 

< kZj M(a)A(a) J J ^-^ (^1 N^tf) J )l 



0{Z b /' 3+s ). 



Proof. We first note that 

(E c (H)+c 2 )(E c (t') + c 2 ) 



(59) 



1 



< 



3E C (£)E C (Z') - c 2 (E c (£) + E e tf) + c 2 ) 



N c {£)Nc{£) 



Then, noting that E c (£) — c 2 < c|£|, wc obtain 



(60) 



(E c (Z)+c 2 )(E c (0 + c 2 ) 



< 



3c 2 |^||ri + 2c 3 (|^| + |C'|) 



< 



2c 4 



Using this last equation, we estimate 

(61) \(F a ,(y--i Pl )F a )\ 



< k I d£d£ 



(c-^HO + ^dfl + lO)) 



/ ~ P) 1 ~ P>> I, :2|,.,l^ , J 



: fc c - 2 i?- 3 / d£ / dg' ^ff ^ffi 9 ^ (II + Rp\\£' + -RpI + cR\£ + iipl) 



< kc~ 2 R- 3 / d£ / d£ 



| C _ € 'ia 
,15(0115(01 



^ . € 7p(l€H0 + «IpII€I + I^pI 2 + cR\Z\ + cR 2 \p\ 

< kc- 2 R- 3 {\ + R\p\ + R 2 \p\ 2 +cR + cR 2 \p\) 
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(see @J Formula (31)]). Thus 

(62) Z\ tr[(| ■ r 1 - tpj-n] <Z\J^ aSl(a)A(a)(F a , (| ■ ^ - ^)F a )\ 

<kZ ( dn(a)A{a) C - 2 R- 3 (l + R\p\+ R 2 \p\ 2 + cR + cR 2 \p\) 



< Z + R / dq Z^VxiZ^q) + R 2 / dq Z^V^Z 1 ' 3 ^ 



5/2 



+ cR + cR 2 J dq yZ^^Vx {Z 1/3 q) 

< k(Z 3S + z 25 + 2 / 3 + Z s+A / 3 + Z 25 + Z s +^ 3 ) 

which yields the desired estimate. 



□ 



3.3.3. The Electron- Electron Interaction. We will roll back the treatment of the 
electron-electron interaction to the treatment of nucleus-electron interaction. 

Lemma 9. For our choice 0/71 and S £ (1/3,2/3) we have 

(63) B(An.. An.) " D (P^Pn) = 0(Z 5 / 3+s ), 

where p 7l is the density 0/71 and p 7llI , is tte density 0/71$. 

Proof. We have 
(64) 

i^[(p 71 + An*) * 1 • r'mi < 2-v\- 5 / 2 \\ P , 1 + P71# H^r 2 = a-vv-v^^i-^. 

Now (see [H Lemma 3]), 

I- (An*' An*) ~ D (Pn>Pn)\ = l-°(An* "An, An* +P7JI 
1 



< 



(AyiW - P 7 i*( x ))[(Ayi + AnJ * I ' I 1 ](x)dx 
< l -J^n{a)A{a) 



x // d£d£'|.F[(p 71 +p 71< 



](€-OI^>OI^(0ll^(O|d€d4' 



<2-3/\s/*Z J dn(a)A(a) J J d£d£'||£ - £r 2 ^(£>O|4(0ll4(O|d£d£' 



where 



(E c (ti)+c 2 )(E c (Z') + C 2 ) 



- 1 



c 2 \m'\ 



and where we used (|64|) in the last step. Eventually, Lemmata [7] and [8] yield the 
desired result. □ 

3.3.4. The Total Energy. Gathering our above estimates allows us to reduce the 
problem to the non-relativistic result of Lieb [21] 

Theorem 3. There exist a constant k such that for all Z > 1 we have 
E(Z/k, N, Z) < E TF (N/Z, 1)Z 7/3 + kZ 20/9 . 
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Proof. Following Lieb [2lJ Section V.A.I] with the remainder terms given there 
(putting R = Z~ s as in our estimate), using the remainder terms obtained in 
Lemmata [7] through [51 and applying (|5"Tj) and (jS"4")) we get 

(65) E(c, N, Z) < ££ F ( 7 ) < £« F ( 7l ) < E TF (N, Z) + 0{Z 1+2S + + zi +s ) 

(see [4j Formula (35)]) which is optimized for 8 — 5/9 giving the claimed result. □ 

3.4. Lower Bound. The lower bound is - contrary to the usual folklore - easy. 
As we will see, it is a corollary of S0rensen's [55] result for the Chandrasekhar 
operator and an estimate on the potential generated by the exchange hole |26j . 
The exchange hole of a density a at a point x g I 3 is defined as the ball Br^ ( x ) (x) 
of radius R a (x) centered at x where R a (x) is the smallest radius R fulfilling 

The hole potential L a of a is defined through 

(67) La (x) := / T^Vy- 

3.4.1. L°° -Bound on the Exchange Bole Potential. We begin with the following 
remark: the Thomas-Fermi potential Vz ■= Z/\ ■ \ — pz * \ ■ | _1 can be written as 
([36"1) . This equation yields immediately the upper bound 

(68) p z (x) < (Z/ 7TF ) 3 / 2 |x|- 3 / 2 . 

This bound allows us to prove the following L°°-bounds on potentials of exchange 
holes. 

Lemma 10. 

\\L pz \\ oo = 0(Z). 

Proof. The function 

/ : R + -» R 

(69) 



t^Vtf lyj-^ly -h (o,o,i)|- 3 / 2 dy 

J|y|<l/t 



'|y|<V* 

is obviously continuous on (0, oo). Moreover, f(t) tends to a positive constant for 
t and to for t -> oo. Thus, ||/||oo < oo. 
This allows us to obtain the desired estimate: 

(70) L Pz (x) < Ai(x)+A 2 (x) 
(see [H Formula (46)]) where 

(71) A lW :=/ £^±Z) dy <(Af7 

7| a |<i/z |y| \7tf / J\ y \<i/z \y\\y + w /2 

= (Z/ 7TF ) 3 / 2 Z- 1 /2/(|x|Z) < ||/||oo7tf /2 ^ 

and 

(72) 4 2 (x) := / Pz( ^t y) dy < Z / Pz(x + y)dy < f . 

•/^<|y|<« PZ (x) ^i<|y|<fl PZ (x) * 

These two estimates prove the claim. □ 
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Lemma [TOl allows us to estimate the N electron operator B c .n.z by the canonical 
one particle Brown-Ravenhall operator whose nuclear charge is screened by the 
the Thomas-Fermi potential. However, since we would like - because of mere 
convenience - to take advantage of S0rensen's result 28 , we derive an estimate 
on L ps (where ps :— pz * g z - S ), i.e., the exchange hole potential of the density 
occurring in S0rensen's proof. 

Lemma 11. 

\\L ps \\ oa = 0(Z). 

Proof. We proceed analogously to the proof of Lemma [101 

T i / W(x + y) f Ps {x + y) 
(73) L P5 (x)< / — dy+ / — dy 

J\y\<l/Z \y\ Jl/Z<\y\<R PS (x) |Y| 

< / dz ff |_,(z) / Pz(x-z + y) dy + z J ps{x + y)dy 

J J\y\<i/z |y| J \y\<R PS (x) 



< I dzg 2 z s (z)A 1 (x - z) + - < kZ 



Z 
~2 

(see [H Formula (49)]) where we used the definition of the radius of the exchange 
hole from the second line to the third line, the definition of A\ in the next step, 
and in the last step the L°°-estimate ([7T|) on A\. □ 

3.4.2. Lower Bound. 

Theorem 4. 

liminf[E(c, N,Z) - E TF (N 7 Z)]Z- 7/3 > 0. 

Z— >oo 

Proof. Pick 6 > and set ps := pz * 9%-s- Then the exchange hole correlation 
bound [26l Equation (14)] implies the following pointwise estimate 

1 N 

( ?4 ) E |x - x I ^E^*!-!" 1 ^)- L ps (x»)]-D(p 5 ,p 5 ). 

l<fi<iy<N ' M y=l 

Because of the spherical symmetry of g we can use Newton's theorem [27] and 
replace ps by pz in the third summand of the right hand side of (|74|) . Then, by 
Lemma ITT| we get that for all normalized € Qn 

(75) £($) > tr[A+(|D | - c 2 - V S )A+}. - kNZ - D(p z , P z) 

where, for t € M, we set [t]_ := min{t, 0} and Vs — Z/\ ■ \ — ps * \ ■ \ . 

To count the number of spin states per electron correctly, i.e., two instead of the 
apparent four, we use an observation by Lieb et al. [23] Appendix B]: Note that 



(76) A_ = U~ 1 A + U, where U 

Indeed, we have 



1 

-1 



A- = - 1 - , A + = - 1 + — - 

and 



RELATIVISTIC GROUND-STATE ENERGY 



15 



We set X := (|D | - c 2 - Va(x))I 2 , and write 
X 



tr 



X 

Thus 

(77) 2tr(A+ 



A, 



> tr A 











A+ = tr A+ 


















))" 




_ 


(» 




= tr ( 


M 



tr A+J7 



X_ 

X_ 



E7 



tr A 



AT_ 
X_ 



X- 








tr A 



X- 




Since |Dq| = E c (p) and X is a 2 by 2 matrix, we obtain 



X 
X 







A, 



= 2tr(X_). 



(78) tr[A + (| J D |-c 2 -n-(x))A+]_ = tr A H 

> tr (V ( x ° Yj= tr(X_) = 2 tr[£ c (p) - c 2 - Vfc(x)] 

Then 

(79) E(Z/k,N,Z) > 2 tr[£? c (p) - 



V«(x)]_ - D{pz,pz) ~kNZ 



(see [U Formula (42)]) where the last trace is spinless. This connects to S0rensen's 
Equation (3.2) from [35]. It is a fundamental result of [31] that E(Z/n,N,Z) = 
E(Z/k, N c (Z), Z), for any N > N C (Z), where N c is the number of negative particles 
that can be bound to an atom of nuclear charge Z. Considering N c < 2Z + 1 (see 
p~9l Formula (1.2)]), when Z — > oo, we can get 

(80) E(Z/k, N, Z) > 2 tr[£ c (p) - c 2 - Vi(x)]_ - D(p z ,p z ) - k O(Z) Z. 



This result then follows using his lower bound. 



□ 



4. Case II: N/Z = A (Constant) 

The case A > 1 has been already solved in Section 3; A = 1 is solved in so it 
only remains to consider A < 1. 

4.1. Coherent States. This Section is analogous to 13.11 

4.2. Upper Bound. This Division is similar to 13.21 and 13.31 Analogously to For- 
mula ([28]) . we introduce the set of density matrices 

(81) S dN := { 7 G e\t)) |£ c ( P h € e 1 ^), < 7 < 1, tr 7 = N = XZ, } 

where 6 1 (f)) denotes the trace class operators on f); A is a number independent of 
c, N, R, or Z. Defining £^ F as Formula (|29p . and as the same as above, for all 
7 G Sqn, we have Formula (|3"u) . 
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4.2.1. Kinetic Energy. According to [25], we can write p^p Z ^ simply as ptf in the 
following. By concavity, and then analogously to Formula (|53[) . we know 

(82) tr[(£ c (p) - c 2 ) 7 ] <tr(-±A 7 ). 
Similarly to Formula (|54|) . we get 

(83) tr[(E c (p) - c 2 ) 7 ] < ~ 7 tf J PtfWcIx + (tr 7 )iT 2 || V. 9 || 2 . 

4.2.2. External Potential. As the same as Part 13.3.121 we can obtain following two 
lemmata by the analogues of Lemmata [7] and [5] 

Lemma 12. For our choice of 7 = J r dil(a)A(a)\F a )(F a \ and S £ (1/3,2/3) we 

have 

(84) 

< Ztr(^ 27 ) < kZ j^Q(a)A(a) J f ^'^M^ME^A = (Z^). 

Lemma 13. For our choice of 7 and 5 € (1/3,2/3) we have 

(85) IZtrKH- 1 -^]! 

= 0(Z 5 / 3+s ). 

4.2.3. The Electron- Electron Interaction. Similarly to Lemma|Hl we § e t Lemma [141 
below. 

Lemma 14. For our choice of 7 and 5 € (1/3,2/3) we have 

(86) D(Pt,Pt) - D{p^ Pl ) = 0(Z 5 / 3+s ), 
where p 7 is the density of 7 and p 7s is the density of 7$ . 

4.2.4. The Total Energy. By the analogue of Theorem G2 we obtain Theorem [5j 
Theorem 5. There exist a constant k such that we have for all Z > 1 

E(Z/k, N, Z) < E TF (X, l)Z 7/3 + kZ 20/9 . 
4.3. Lower Bound. This Part is similar to 13.41 

4.3.1. L°° -Bound on the Exchange Hole Potential. 
Lemma 15. 

\\L PT¥ \\ 00 = 0{Z). 

Similarly to Lemma fTTl since we would like to take advantage of S0rensen's result 
|28j . we derive an estimate on L ps (where ps := ptf * 9^-^)- 

Lemma 16. 

\\L ps \\ oo = 0{Z). 

4.3.2. Lower Bound. 
Theorem 6. 

liminf[£(c, N, Z) - E TF (N, Z)]Z- J/3 > 0. 

Z—too 
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Appendix A. The Proof of B C: z = E c {p) — Ztpi — Zip 2 

Lemma 17. The one-particle Brown-Ravenhall operator £> 7 for an electron the 
external electric potential of a point nucleus acting on Pauli spinors is 

B c ,z = E c (p) - Zipx - Zip 2 . 
Proof. For any ip £ Sj, we can get 

(87) £{i>) = (il>,B c , z il>) = ^,D c ^) 

- ( Ec(t>)+c\ ci><r A ( <?-Z\-\~ X ™ V )( E ^SfA 

\ JV c (p) N c (p) U ) { ccr . p -c 2 -Z\-\- 1 \ J^L, 



(u, ($ic 2 $i + $iccr • p$2 + <J> 2 ccr • p <f>i $ 2 c 2 $2) - Z(u, <£i«) - Z(u, tp 2 u) 

= (u, E c (p)u) - Z(u, tpiu) - Z(u, (fi2u). 

□ 



Appendix B. Checking of 71 and 7 

Lemma 18. 

< 71 < 1 and < 7 < 1. 
Proof. For any normalized u <E f), using Parseval's equality, we get 

(88) (u, 7l «) 

= (2ir)~ 3 J u(x) J J /(x — q)/(x'— q)A(a) exp(— ip-x) exp(ip-x')u(x')dpdqdx'dx 
= y (2n)- 3/2 J cxp(-ip • x)u(x)/(x - q)dxA(a) 
x (27r)^ 3 / 2 y cxp(-ip • x')u(x')/(x' - q)dx'dpdq 
= y (2vr)" 3/2 y exp(-ip • x)u(x)/(x - q)dx 

< y (2tt)" 3/2 y cxp(-ip • x)u(x)/(x - q)d. 

= yy | U (x)| 2 |/(x-q)| 2 dxdq 
= y Kx)| 2 [y |/(x-q)| 2 dq 



A(a)dpdq 

2 

dpdq 



dx = llull = 1, 
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and following 
(89) (u, ju) = (u, fit*) + (u, 7 2 u) + (u, f 2 u) 

K+N-Z 

< 



J | U (x)i 2 dx+ <r J | U ( X )i 2 dx 

\x\<R-R k=K+l xesupp{0 s } 

K+N 

+ e y i«(x)i 2 dx 

fc=if+JV-Z+l xesupp {0 s } 

< Ilu|| a = l. 



Appendix C. Checking of tr7i and tr7 



Lemma 19. 



tr 71 < N and tr 7 = JV 
Proof. Using formula (1361) . we can get 

2 



(90) tr 7 i= J m{a)A{a) = (2tt)- 3 / ^ dpdq 

' 1 p 2 /2-Vz(q)<0 T=1 

= 2(2^)- 3 ^(2 7TF ) 3 / 2 y p z (q)dq = Z < N, 



and 



(91) tr7= J an(a)A(a) + (N - Z) + e R Z = N. 



□ 



□ 
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